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A system of metastable plus unstable states is discussed. The mass matrix 
governing the time development of the system is supposed to vary slowly with time. 
The adiabatic limit for this case is studied and it is shown that only the metastable 
states obtain the analogs of the dynamical and geometrical phase factors familiar 
from stable states. Abelian and non-abelian geometric phase factors for metastable 
states are defined. 
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1. Introduction 

The investigations reported in this paper are in connection with an effort to measure 
atomic parity- (P-) violating effects with an atomic beam interferometer described in 
Parity violation in atoms through neutral current exchange was already studied in 2fl. The 

n n n n 

seminal paper in this field is |3|], for reviews see |4|, [5|, |6j. Measurements of P- violating 
effects in heavy many electron atoms like Cs and Tl have already reached a high precision 
5, 7, f|. For atoms with a single electron like hydrogen or deuterium, however, experimental 
measurements have, so far, not been successful. In an atomic beam apparatus an atom can 
be subjected to external electric and magnetic fields. Under suitable conditions the motion 
of the atom can be described in the adiabatic limit. As is well known from Berry's work 
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jsl the atomic wave function can then acquire a geometric phase factor. For a collection of 
important papers on geometric phases see [10]. We shall study how such phase factors occur 
for metastable states and classify the phases in parity conserving (PC) and parity violating 
(PV) ones. The final aim is to see if PV geometric phases are large enough to be measurable 
with an atomic beam apparatus. But this will be dealt with in future work. 

Our paper is organised as follows. In section [2] we discuss the description of unstable 
atoms in slowly varying external electric and magnetic fields. We discuss the adiabatic limit 
for stable states and a set of unstable states all having the same decay rate in section 3. 
In section 4 we deal with a simple example of a two state system where one member is 
longer lived than the other one. The more general case of a number of metastable states 
with equal decay rates and a number of unstable states with larger decay rates is treated in 
section 5. Our conclusions are presented in section 6. Two appendices present the detailed 
mathematical arguments leading to the results of sections 4 and 5. In the companion paper 



II [llj we study the metastable 2S states of hydrogen and deuterium and identify the PC 
and PV phases occurring there. In the following we shall refer to tables of II by table II. 1 
etc., to equations as (II. 1.1) etc. 

We use units h = c = 1 if other units are not explicitly indicated. 



2. Metastable states, Generalities 

Let us consider an atomic system with N stable or unstable states. We assume that the 
atom is at rest and subjected to slowly varying external electric and magnetic fields. The 
effective Schrodinger equation for the system is then, in the Wigner-Weisskopf approxima- 
tion, 

i^\t) = MW) , (2.i) 

where \t) is the state vector of the undecayed states and ^(t) the mass (or complex en- 
ergy) matrix which is, in general, non hermitian. For a discussion of the complete Wigner- 



Weisskopf solution see for instance 121 ] . 



To give an example which we shall study in detail in II let us consider the states with 
principal quantum number n = 2 in ordinary hydrogen. As basis states we have here the 
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25" and 2P states. In total these are iV = 16 states which we number as shown in table 
II JU In table IIj5]we also give the mass matrix for the case of zero external fields. With 
external fields the mass matrix is 

J£{t} =JK q -D- £(t) - n ■ B(t) ■ (2.2) 

Here £{t) is the electric, 3{t) the magnetic field strength vector. These are supposed to 
vary only slowly with t. Furthermore .D is the matrix of the electric and /x of the magnetic 
dipole operators in the space of n = 2 states, see tables IIJ6] and 1113, respectively. 

We return to the general case (12. ip . We want to study the adiabatic limit, that is, a given 
change of with t is made over a longer and longer time. Following [l3] we implement 

this by introducing a reduced time r, setting 

t = — t (2.3) 

where r is some fixed time. We consider the system over the interval 

< r < t (2.4) 

in the limit of larger and larger T. The mass matrix in (12.11) is supposed to be a function of 
r only 

^L{t) = Mr) ■ (2.5) 

We get then from (O) with \t) = \T;r) 

\^\T-t) = -Mj)\T-t) . (2.6) 
or r 

Here and in the following we indicate explicitly the dependence of the quantities on r and 
T, respectively. 

The mass matrix ^(t) ( 12.51) can be decomposed into the hermitian energy and decay 
matrices 

Mr) = E{t) - i£(r) , 

E}(t) = E(t), ( 2 - 7 ) 

rt( T ) = r(r) . 
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We suppose that £(r) is a positive-semidefinite matrix 

£(r) = \{Jttf) - J£{t)) > . (2.8) 

We suppose, furthermore, that for all r ^#(r) can be diagonalised. (This is, for instance, 
guaranteed if j&{t) has N different eigenvalues.) Then we have for all r complete sets of 
right and left eigenvectors satisfying 

^#(r)|a, r) = E(a, r)\a, r) , (2.9) 

(c^ri^(r) = (5rr|S(a,r) , (2.10) 
(« = l,...,iV) . 

Here E(a, r) are the — in general complex — eigenvalues of jM{t\ 

E(a,T)=E R (a,r)- 1 -T(a,r), (2.11) 

with Eft(a, t) the real part of the energy and T(a, r) the decay rate of the state \a, r). We 
choose the normalisation of the eigenvectors such that 

(atf\P,T) = Sap , (2.12) 
(a, r\a, r) = 1 (no summation over a) . (2-13) 

Everything is supposed to be continuous in r and we shall use this in the following in an 
essential way. 

From (US} follows 

r(a,r)>0, (a = l,...,N). (2.14) 

The proof is simple. From (I2.8l) - fl2.13l) we obtain 

(a,r\T(r)\a,T) = (a,r\i L4{t) -^{t) \ \a,r) = T(a,r) > . (2.15) 

Note that the reverse is in general not true, that is, (12.81) does not follow from (12.141) . 
Now we expand the state vector \T\ r) in terms of the eigenstates \a, r): 

N 



\T;r) = J2MT;r)\a,r) . (2.16) 



a=l 



From ( 12. 6p we get then 



i^-MT; r) = -E{a, r)^ Q (T; r) - ^(^|i|-|/3, r)^(T; r) . (2.17) 



The task is to discuss the solution of (I2.17P for large T 
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3. Stable states and states with equal decay rates 

Here we discuss briefly the cases that all states are stable or that they are unstable but 
have equal decay rates. In detail we suppose the following 

r(a,r)-r(/5,r) = (3.1) 

for all a, (3 G {1, . . . ,N} and all r G [0,r ]. This includes the case of stable states where 
r(a,r) = 0. Furthermore we suppose 

\E(a,r) - E((3,t)\>c> (3.2) 

for all a ^ [3 and all r G [0,r ], where c is a positive constant. With (13.11) and (13. 2p the 
discussion of the adiabatic limit of the solution of (I2.17P can be done exactly as in 13( and 
one finds as solution of (12.171) for large T 



V»«(T;0)+O(i 



(3.3) 



r t 

ip a (T; t) = exp -i— ip a (r) + i-f aa (r) 
L r o 

(a = l,...,N) . 

Here T<^ q (t)/to is the dynamical and ^ aa {j) is the geometric (Berry-)phase, 

<P«(t) = [ dr' E(a,r') , (3.4) 
Jo 

r -- d 

l aa (r) = / dr (a,r'\i— \a,r) . (3.5) 



dr' 

From (13. 3p we see that for large T the following holds. If we start with an eigenstate of 
^#(0) for r = the system will be in the corresponding eigenstate of ^# (rn) for r = ro. 
For decaying states (satisfying (13.11) ) both, the dynamical and the geometric phase factors 
will have real and imaginary parts. That is, also the effective decay rates of the states a 
will obtain a geometric contribution. The real, but not the imaginary part of the geometric 
phase is "gauge" dependent. If we change the definition of the states \a,r) by 

\a,r) — ► \a,r)' = e iVa( - T) \a,T) , 



(3.6) 



where r) a (r) must be real in order to respect (12.131) we get 



9 .. ,. d , dr] a (r) 

{a,T\i— \a,T) — ► (a,r\i—\a,T) . (3.7) 
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The change of the geometric phase (I3.5P induced by the transformation (13. 6p of the basis 
states is, therefore, 

l' aa ( T ) = laa{r) - r) a (r) + r) a (0) . (3.8) 



4. A two state system of decaying states 

Before going to the general case of N decaying states we discuss as an example a system 
with two unstable states (N = 2). Let the state 1 be the longer-lived one. In detail we 
suppose 



r(2,r)-r(i,r) > Ar min >o 



(4.1) 



for all r G [0,r ], where Ar min is a fixed constant. The real parts of the energies, Er(1,t) 
and Er(2,t), see ( 12. lip , can be arbitrary. From (12. 17ft we get for N = 2 the coupled 
equations 



T 

TO 

T 

TO 



-E(1,t) - on(r) 



E(2,r)-a 22 {r) 



^i(T; r) -Oia(r)V^(T;r) 
ip 2 (T\ r) -a 2 i(r)^i(T;r) 



(4.2) 



Here and in the following we use the definitions 



d 



a a p{r) = (a,r\i— \(3,r) 



7o/j(t) 

(pair) 



dr' a a p{r') , 
dr' E(a, r') 



dr' 



E R (a,r')--T(a,T r ) 



Pa{r) 



dr' T(a, r') 



(4.3) 
(4.4) 

(4.5) 

(4.6) 



where a, (3 G {1, 2}. (In later chapters we use these same definitions for a, (3 G {1, . . . , iV}) 
We have 



Imyj Q (r) 



1 



Pa(r). 



(4.7) 



With the initial conditions ip a (T;0) we can transform (14.21) into integral equations 



ipi(T] r) = exp 



T 

-i— (pi{T) + i7u(r) 

TO 



x \ ipx(T;0) + / dr' exp 



T 



-i— y? 2 (r) + i72 2 (r) 

7b 



ip 2 (T;r) = exp 

x <U 2 (T;0) + / dr' exp 



T 

i— <£>i(r ) - i7n(r' 

To 



i— <^ 2 (r ) - i722(r ) 

To 



ia 12 (r')^ 2 (T; r') 



(4.8) 



ia 2 i(r')^i(T; r') 



(4.9) 



Inserting (14.91) in (14.81) we get an equation for ipi(T; r) alone. To write this in a transparent 
form we introduce functions Xi(T; r) and Xi (T', r) by 



^i(T;r) = exp 



T 

-\—<Pi{t) + i7n(r) 
7) 



7i (T; r) 



xS 0) (T; r) = Vi(T; 0) + Z 12 (T; r)^ 2 (T; 0) 



(4.10) 
(4.11) 



where 

Ii2(T;t) = i / dri exp 



T 
i — 

r 



{(pi{n) - <P2(ti)) - 1711(71) + 1722(71) 



a 12 (n) . (4.12) 



Furthermore we define an integral operator L for continuous functions r — > £(r) (r G [0, r ]) 
as follows 

T 



MO 



dri exp 



x / dr 2 exp 
'0 



(<£i(ti) - ^(n)) - i7n( r i) + 1722(7-1) 
T 

-i— (</>i(t 2 ) - ^2(7-2)) + 1711(7-2) - 1722(7-2) 

7"0 



ai2(Tl) 



021(7-2) ((t 2 ) . 



(4.13) 



We get then from (Oil and (Oil 



7 1 (T;r) = 7j 0) (T;r) + (L Xl )(T;r) 



(4.14) 



From (I4.14p we obtain 



((l-L)7 1 )(T;r) = X f ) (T;r) 



(4.15) 



with the solution 

00 

Xi(T; r) = ((1 - L)^x! 0) ) (T; r) = X f \T, r) + £ (W° } ) (T; 
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(4.16) 



n=l 
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With fl4XL]) this gives 

oo 

Xi(T; r) = ipi(T; 0) + Z ia (T; r)^ 2 (T; 0) + ]T (l» X <°>) (T; t) . (4.17) 



n=l 



Now we show that all terms on the r.h.s. of (14.171) except ipi(T; 0) are of order 1/T for 
large T. We note first that from (14.51) and (14.61) we have 



exp 



T 

±i—<p a (r) 
To 



exp 



1 T 
±-—p a (r) 

2 T 



(a = 1,2) 



(4.18) 



We get now from (JQ2l) and (JO 

|ii 2 (T;r)| < / dri exp 



IT 
2^ 



(p 2 (ri) -pi(ri)) 



exp ( - i7u(ri) + i72 2 (n)) 1 012(71) I 



< / dri exp 
'0 



1 T 



2 T 

x(r(2,r 1 )-r(l l ri))Ari 



(p 2 (ri) -pi(ri)) 



(4.19) 



exp ( - i7n(ri) + i7 22 (ri)j |ai 2 (ri)| 



Above, after (12.131) . we have explicitly supposed that all functions like Oi 2 (r), 7n(r) etc. are 
continuous for r G [0, r ]. Thus they are bounded there. We set 



C12 = max <5 t 

ne[0,r ] 



exp ( - i7n(ri) + i7 22 (ri)j |ai 2 (ri)| 



(4.20) 



From (I4.19P we get then with (j4T6|) 



Vi2(T] t)\ < ci 2 I dn exp 
Jo 



(P2(n) -Pi(n)) 



= C12 / dri exp 
2 / 

= C12- II - exp 



T 
'2T0 
T 

■— (P2(ri) -pi(ri)) 
zr 

T 

7T- (p2(r) -pi(r)) 
zr 



[T(2,ri)-r(l,ra)] (r Ar mi 
^r(p2(n) - Pi(r 1 ))(r Ar E 

J (Ar^)- 1 



(4.21) 



Here we use that due to (14.11) and (14. 61) we have for all r 6 [0, r ] 



p 2 (r)-pi(r) >rAr min >0 . 



(4.22) 
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In appendix |A] we show that for large T 



CO 
71=1 



< Ci(TAr min ) _1 |^(T; 0)1 + C 2 (T Ar r 



\MT;0)\ 



(4.23) 



where C\p are positive constants. We get, therefore, from (14.171) with (14.211) and (14.231) for 
large T 



Xi(T;r) = ^i(T;0) 



+ C9(i)r 2 C/':0). 



(4.24) 



Inserting this in (14.101) we find that for large T the solution of the original amplitude for the 
longer-lived state is 



r) = exp 

x (v>i(T;0) 



T 

i— <Pi(r) + i7n(T) 
To 



(4.25) 



For the amplitude of the shorter-lived state we obtain, in an analogous manner, from ( 14. 9 h 
and (14.251) (see appendix [AT) 



ip 2 {T; t) = exp 
+ exp 



T 

-i— <Pi(r) + i7n(r) 
T 

-i— <£ 2 (r) + i722(r) 

TO 



1 

y2 



(4.26) 



The interpretation of the results (14.251) and (14.261) is clear. The amplitude for the longer 
lived state, a = 1, shows the usual dynamical and geometrical phase factors, Ty9 1 (r)/r 
and 711(1"), respectively. Both these factors are in general complex for decaying states. The 
shorter lived state, a = 2, gets in general some feeding from the longer lived one. Thus the 
leading term of its amplitude (the first on the r.h.s. of (14. 261) ) has the same exponential 
factor as the amplitude for a = 1. But there is a suppression by at least a factor of order 
1/T. For r > and large T the second term on the r.h.s of (14.261) vanishes exponentially 
relative to the first term, see (14.180 and (14.221) . 

In particular, if we start at r = with ipi{T] 0) 7^ and ip2{T; 0) = 0, that is only the 
longer lived state is populated, the system stays in the state a = 1 up to corrections of 
relative order 1/T. The amplitude evolves with the appropriate dynamical and geometrical 
phase factors. If, however, we start at r = with ip\(T; 0) = and ip2{T; 0) 7^ 0, that is only 
the shorter lived state is populated, it will in general populate to a small extent the longer 
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lived one. The latter will then feed back on the former and "impose" its decay properties 
on the system. But all this gives only amplitudes which are suppressed for large T. To 
summarise: the adiabatic theorem in the usual form can only be used for the amplitude of 
the longer-lived state. 



5. Several metastable and short lived states 



Here we generalise the results of section H] to a system having several metastable and 
several short lived states. In detail we suppose that the system has N states where the first 
M ones (1 < M < N) are metastable. We suppose equal decay rates for the metastable 



states, that is 



r(a,r)=r(l,r) 



(5.1) 



for a G {1, . . . , M} and r G [0, ro]. The remaining N — M states should have substantially 
larger decay rates for all r: We suppose 



T(P, t) - T(a, t) > Ar min > 



(5.2) 



for all a, f3 with 1 < a < M and M + 1 < (3 < N and all r G [0, r ]. Here Ar min is a fixed 
constant. 

The problem is again to solve (I2.17P for large T. As in section H] it is convenient to 
transform (12.171) into an integral equation. We use the same notation as in (14. 3ft - (14. 71) but 
now for 1 < a, j3 < N. With the initial conditions ip a (T; 0) we obtain from (12.171) 



T 

i — iPa[T) + ij a a{T) 
TO 



ip a (T; t) = exp 

x |vo(T;0) + J dr' exp 
(1 < a, (3 < N) . 



i—ip a (T) - ij aa (T') 



In the following we use matrix notation. We set for the metastable states 



ip a (T; r) = exp 
/ 

x(T;r) = 



T 

-i—ip a (r) + i7 QQ (r) 



X«(T;t) , (l<a<M) 



Xi{T; r 



\ 



\Xm(T; t)J 

For the short lived states we set with pi(r) from (14.61) 



T 



^p(T; t) = exp 



^(T;r) , (M + 1<P<N) 



We have for r = 



\ 6v(T;r) J 
( 



\ 



= X (0 \T;t), 



^M + l(T;0)^ 



e(T;0) 



Furthermore we define 



£(°)(T;r) 



V <MT;0) J 



f(°)(T;r)=exp 



T T 
-i— <f a (r) + — pi(r) + i7 aQ (r) 
To 2r 



^a(T;0) , 



(M + 1 < a < N) . 
The integral equations (I5.3P can then be written in matrix form as follows 

X = X i0) +L^ X + L^, 



12 



Here the L^'fi are integral operators (generalising ( 14. 12ft and ( 14.13P ) which are given explic- 
itly in appendix El 

The formal solution of (15 . 1 2[) and (15.131) is easily written down. We get from (15 . 13[) 



(l-L^)£ = ^ + L^x, 



(5.14) 
(5.15) 



Inserting this in ( 15.12ft gives 

|l _ _ L (l,2) _ L (2,2))-l L (2,l)| x = x (0) + L (l,2) ^ _ L (2,2))-l |(0) ? (5 Jg) 



x = {i - - ^ (1>2) (i - L^y 1 l&v} 1 { 



X 



(o) 



1_ L (2,2))-1|(0)| (517) 



This is the exact solution. In appendix [B] we show that for large T we get from (I5.17P the 
result 



X 



(5.18) 



or written in another form 



x(T,T) = U(T- 1 r) X ^ + 



T 



) llx (O) ll+0(^) ll^ (0) ll , (5-19) 



see (1BT82]) and (lR83|) . Here the norms ||x (0) l|, ll£ (0) ll and the M x M matrix £/ are defined 
in (|R3|) . (lR4j) and (lBTl6|) . respectively. 

The solution for £(T; r) is given in (15351) . The term (l-L^ 2 ))- 1 ^^ 1 )^ can be estimated 
using (1B.55|) and we find 



(5.20) 



With (15TI9|) . (|B~T6|) and (iRlB we get from (CTj) 

||(1 - L P, a ))-l L (2,D x || < O || X (0)|| + ^ || e (0)|| . (5 . 21) 



The term (1 — L^ 2 ' 2 ^) 1 ^°> in (15. 15j) can be estimated according to (1B.71I) 



T 

2^" 



AT ■ t 



lle 0) || ■ 



(5.22) 
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Putting everything together we find 



\X 



(0)1 



O 



1 

2^2 



lle (0) ll + 



exp 



T 



-Ar ■ r 



(0)1 



(5.23) 



Now we can go back to the original amplitudes tp a (T; r) of (I2.17P and (15. 3p . We find from 
(15. 4p . (15. 8p and (15.191) for the amplitudes of the metastable states 



ip a (T; t) = exp 



T 

-i—cp a (r) 

TO 



X 



./3=1 



(5.24) 



for a G {1, 



u5T(r;r; 



(5.25) 



, M }. Here the matrix of the geometric phase factors is given by 

U geom (T;r) : 

U^ m (T; t) = exp[i-f aa (r)]U af3 (T; r) 

with U (T; r) given in (IB. 16[) . For the amplitudes of the short lived states we find from (15.61) . 
(153|) and (l5"^3l) 



T 



exp 
+ exp 
P G {M + 1 



T 



\X 



10111 -o(^)ne (0)l 



-— (px^+Ar^nr) 
,iV}. 



(0)1 



(5.26) 



The interpretation of (I5.24p and (I5.26P is analogous to the one given for the two-state 
system in section 4. Suppose that for t = only the metastable states are populated, 
that is, we have an d C = 0. Then the amplitudes for the metastable states 

evolve according to (15.241) . obtaining a dynamical phase factor but also a — in general non 
abelian — geometric phase factor given by the matrix U geom (T\ r). Corrections to this are 
suppressed by a factor T _1 for large T. Looking through the formulae of appendix [B] we see 
that these suppressed terms are of order (TAr min ) _1 relative to the leading term. We shall 
analyse the metastable states further in section [61 

On the other hand, suppose that initially only the short lived states are populated, that 
is, we have = an d ^ 0. Then the amplitudes of the short lived states have one 
part showing a fast decay corresponding to the decay rates of the short lived states and 
terms 0{1/T 2 ) showing the decay as for the metastable states. The metastable states get 
in this case only amplitudes suppressed by 0(1/T). 



14 



6. The metastable states 



In this section we study the evolution of the metastable states alone. That is, we assume 
T to be large enough such that all terms 0(1/T) can be neglected. We have then from 

M 

MT; r) = J2 U a /3{T; r) X f , a e {!,..., M}, (6.1) 



0=1 



where 



U{T ] r)={u aP {T-r)) 



U dyn (T;T)U 8eom (T;r) , 



U dyn (T- 



t) = exp 



T 

T 



>afi 



and m com {T; r) is defined in fl5T25|) . This gives with U from (15161) 



U a f3(T; r) = exp 
From (g3D, (USD, flRTBl and flRTTl) we find 



T 

To 



U a p(T;r) . 



(6.2) 
(6.3) 



(6.4) 



U{T;0) = 1 



M 



d_ 



U{T-t) = -i^(T;r)U(T;r) 



where 



T 

^,g(T; r) = — £(a;, r)5 Qj g - a a/3 (r) 
a,/3e{l,...,M}. 



(6.5) 
(6.6) 



(6.7) 



Thus, the amplitudes for the metastable states ( 16.11) evolve according to the effective 
Schrodinger equation 

M 

(6.8) 



i-^MT; t) = J2 ^M T 5 t )MT; t), a G {1, . . . , M} . 

0=1 



with effective mass matrix JV (T; r). 
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Recall that we supposed equal decay rates for the metastable states. If the energy eigen- 
values for these states satisfy (13. 2p we can use the results of section |3j Then the amplitudes 
ip a {T; t) (16.81) decouple in the evolution for T — > oo and we get up to corrections of relative 
order 1/T (see (13~3D ) 



T 



-i—<p a (r) + i7 aQ (r) 
To 



^ a (T; r) = exp 
U(T;t)= (exp 
U^(T;t)= fexp [i laa (r)) 5 aP 



y (o) 



T 

i—<p a (r) +h aa (r) 

TO 



>aj3 



(6.9) 

(6.10) 
(6.11) 

That is, in this case the matrix of geometric phase factors becomes diagonal. In the general 
case, however, JJ geom (T; r) will contain non-diagonal terms. 

Finally, we note that the results of this paper remain the same if the metastable states 
have slightly different decay rates and we consider only times T for which 



[l>,r)-r(/3,r)]T|«l 



(6.12) 



for all a, P e {!,..., M} and all r £ [0, r ]. 



7. Conclusions 

In this article we have analysed the temporal behaviour of a system consisting of several 
metastable and several short lived states. The mass matrix governing the evolution of the 
system was supposed to be slowly varying with time, see (I2.1l) - (l2.5p . The adiabatic theorem, 
adapted to decaying states, was shown to hold for the metastable states, see section [5j The 
evolution of these states (see (16.11) . (16.21) ) is governed by a dynamical phase factor matrix 
U yn (T,r) (16.31) which is diagonal and a geometric phase factor matrix U gcom (T;r) (15.251) 
which, in general, is not diagonal. 

In the accompanying paper II we apply these results to a study of the states of hydrogen 
and deuterium with principal quantum number n = 2 in slowly varying external electric 
and magnetic fields. The 25' states are metastable and we shall define and discuss parity 
conserving and parity violating geometric phase factors for them. 
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Appendix A 



Here we give the details of the estimates (14.231) and (14.261) for the two-state system 
discussed in section HJ 

We start with the proof of (14.231) . The operator L is defined in (I4.13p . For any continuous 
function r — > ((t) on the interval [0,Tq] we have with (14. ip and (14.181) 



|(LC)(r)| < / d n exp 
Jo 



n 

x / dr 2 exp 
'o 



T 

2^(P2(n) -Pi(ri)) 
T 

^r(p 2 (T- 2 ) - Pi(r 2 )) 



|exp [-i7n(r 1 ) + i722(ri)J| 



[r(2,r 2 )-r(i,r 2 )] 



x Ar min l ex P [1711(^2) - i722(r 2 )]| |a 2 i(r 2 )| max |C(r')| • 

T'e[o,T ] 

Analogously to (14.201) we define 

c 2 i = max {ro|exp[i7u(r 2 ) -i7 22 (r 2 )]| |a 2 i(r 2 )|} 
T2e[o,ro] 

where c 2 i is a finite, non-negative constant. We get then with (I4.20p and (1A.2P 



{L(){t)\ < c 12 c 2 it Ar min max |((r')| / dn exp 



Tl 

x / dr 2 exp 
'0 



e[o,r 
T 



2t, 







(p 2 (r 2 ) -Pi(r 2 )) 



dTn 



T 

-7T (P2(ti) -Pi(ti)) 
(p2(t" 2 ) - Pi(r 2 )) 



c 12 c 2 i r 2 Ar m J n max |C(r')| 

T'e[o,r ] 



x — J dn <( 1 - exp 



T 

-tt~ (Piin) -pi{n)) 



< c 12 c 21 ATj m - max \((r')\ 

1 T'e[0,r ] 



Setting 



(A.l) 



(A.2) 



(A.3) 



t = 2ci 2 c 2 i Ar_L 



(A.4) 



we have thus 



|(LC)(r)|<^ max |C(r') 

1 1 1 r'e[0,T ] 



and by straightforward iteration 



r 



\(L n 0(r)\<(^\ max |C(r')| 

1 1 \1 J r'e[0,-ro] 

For large enough T we have certainly T > 2fo and we get 



n=l 



< 



£ 

n=l 



| max |C(t')I 



= 5(1-5) max \C(t')\ 
< ^ max |C(r)| . 



Replacing here £(t) by xS 0) (T; r) we get with fTCTT]) and EI} 



^(L^ 0) )(T;r 



n=l 



. 2f 

< max 

T r'e[0,r ] 



xS 0) (T,r') 



<Y-{\MT',o)\ + o(^)\MTv 



which proves (14.231) . 

For the proof of (14.261) we start from ffl~9j) and (I4.10P to get 



MT;r)=iJ ( i\T ] r) + ^>(T;r) , 



(2), 



where 



r 



exp 
exp 



T 

-i— ^i(r) + i7u(r) 
tq 

T 

-i— (^ 2 (r) - ^i(r)) 



exp [-i7u(r) + 1722(7*)] 



x / d,T\ exp 
Jo 

xXi(T;n) , 



i — 

T 



(^ 2 (ri) - <pi{n)) 



exp [i7n(ri) - 1722(71)] ia 2 i 



exp 



T 

-i— <p 2 (t) + i722(r) 
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Using the same techniques as for (L() (r) we find with c 2 i from (1A.2j) and with (14.241) 



<c 21 (t AT 



mm , max |exp -1711 (r ) + 1722(7 

t'£[0,t ] 



T 



x exp 

x / dri exp 



- — {p2(r) ~ Pi(r)) 



T 



2 



2r 
l + O 



(p2(n) -pi(n)) 
+ 



(P2(ri) -Pi(ri)) 



l^(T;0)| 



(A.13) 



< ^ Ar min c 2i max J exp [-i7n(r') + i7 22 (r') 

1 r'e[0,r ] 



x |^(T;0)| 



O ( ^ ) l^(T; 



Inserting (lA~T0l) - flA~T3l) in ([Oil proves (T4T26D . 

To summarise, we find the following results where for later use in II we write out the 
nominal order of magnitude of the first correction terms 



■0i(T;r) = exp 



T 

-i— <fi(r) + i7n(r) 



^ 2 {T; t) = exp 



x V>i(T;0) 



T 



l + i^L (i) 

TAr ■ 

J. z_\l min 



+ ^^<WMT;0) 

^ min 



(A. 14) 



-1— v?i(r) +i7n(r) 
2 '" 2J 0(l)Vi (T; 0) + J C A 2 l Cl \, 0(1)MT; 0) 



TAr ■ 



+ exp 



T 

-i— <P 2 (t) +i722(r) 
t 



(TAT mi 



(A.15) 



Appendix B 



Here we give the details of the calculations in section 

Inserting (15.41) - ( T577T) in (15.31) we get (15.121) and (15. 131) where we define the operators L^> 
'1 < i,j < 2) as follows. Let 



7--C(l 



\Cm(t) y 



(B.l) 



and 



^m+i(t)^ 
T -> 77(r) = i 

be continuous vector functions for r G [0, r ]. We define the norm for such functions 



= max |C a (r)| , 

l<a<M, 
t€[0,t ] 

= max \r}p{r)\ 

M+l<(3<N, 
re[0,T ] 



The operators L^'fi are defined as 



T 

I — Vat?) ~ i7aa(T / ) 
70 



(^ (1 ' 1) C) a (r) = i /V exp 

./0 

M 

x ^a a/3 (r') (1 - <5 a/3 )exp 



0=1 

(1 < a < M) , 



T 

-i— <Pp(t') + i7/3/3(r) 

TO 



{L^ 2) v)Jr) = i f dr' exp 
Jo 



T , 
i—ip a (r') - haa(r') 

TO 



N 



x ^ a a/3 (r')exp 

P=M+1 

(1 < a < M) , 



T 
2r 



W(0 , 



(L^C) a (r) = exp 

x i / dr' exp 
Jo 



T T 

i— <A*(r) + — pi(r) + i7 aQ (r) 

7"0 ^0 

i—(p a (r') - haa(r') 

TO 



M 



x ^a a/3 (r')exp 

0=1 



-i— (pp(r') + i7/3/3(r) 

T"0 



(M + 1 < a < iV) , 
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v rj = exp 



T T 
-i—<Pa(r) + — pi(r) + i7oa(r) 
to 2r 



T 

i—(f a (r) - ijaa(r 

TO 



x i / dr' exp 
Jo 

N 

x a a ^(r') (1 - 5 af3 ) exp 



(B.8) 



T 



^(t) . 



(M + I < a < N) . 



Our aim is to derive estimates similar to f)4.2ip and flA.3p ~ flA.8p in order to prove 05.180 . 
Let us first recall that we are supposing all functions of r to be continuous for r G [0,r ]. 
Thus we can define the following dimensionless nonnegative and finite constants 

a^fY) (1 - 5 af3 ) exp [-i7aa(r) + il^{r)] 



C\\ = max t 

l<a<M, 
1</3<M, 
tG[0,t ] 

C\2 = max t 

l<a<M, 
M+l<0<N, 
re[0,r ] 



C21 = max r 

Af+l<a<iV, 
1<0<M, 
re[0,T ] 

C 22 = max t 

M+l<a<N, 
M+1</3<N, 
re[0,To] 



a a(3 (r) exp [-i7 aa (r) + iT^C 7 ")] 

a aj3 (r) exp [-i7 QQ (r) + iT/^Cr)] 
a a/3 (r) (1 - (5^) exp [-i7 aa (r) + vypp{r)\ 



(B.9) 
(B.10) 

(B.H) 
(B.12) 



Now we consider the solution for x given in f)5.17p . With simple algebra we can rewrite 
it as follows 



X = (1 - L^))- 1 (l - I)' + (1 - L^y 1 

= (1-L^y 1 (f2 Lk ) {* (0) + L(1,2) (1 -^ (2,2) ) _1 ^ (0) } 



(B.13) 



where 



Z = L (1 ^ (1 - L^ 2 ))" 1 L (2,1) (1 - L^)' 1 . 
We study first the operator (l — L^ 1 ' 1 )) . Let 



(B.14) 



(B.15) 



21 



be a constant vector. We define a M x M matrix function U(T; r) by 

U(T;t)( c = (fi-L^y 1 ^) (r) 

oo 

= E(( i<M) ) r Co)(r) 



r=0 



Cc+ ^Et^TCc (r) 



r=0 



From flB.5j) and (IB. 16|) we find 



d_ 



d_ 



T 



r=0 



A(T;t)U{T;t)( c 



where 



A(T;t) = (A aP {T-T)) , 
A a p(T; t) = exp 



T 

i— y?a( r ) - i7c«*0~) 

?0 



ia a p(r) (1 - 5 a/3 ) 



x exp 
(1 < at,P < M) . 



-i— <Pp(t) + vyppif) 
To 



We have 



U(T;0) = t M , 
det U{T; 0) = 1 , 

Tr r) = , for all r G [0, r ]. 



We also have the relation 



9 

— det U(T; r) = Tr ACT; r) det C/(T; r) 



following from flB~T7l) . We find, therefore, with flR20|) and flB~T9l) 



det *7(T; r) = det t>(T; 0) = 1 



(B.16) 



(B.17) 



(B.18) 



(B.19) 
(B.20) 

(B.21) 

(B.22) 



for all t G [0,ro]. Thus U(T;t) is never singular and has an inverse U 1 (T;r) for all 
t G [0, To]. For the inverse we have 





Or 



U-\T- t) = -U- 1 (T;t)A(T;t) 



(B.23) 
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The matrix elements A a ^(T; r) (IB. 181) are bounded as we see from (1B.9j) . (I4.18P and (15.11) : 



\A a p(T;T)\ < 



T 



(B.24) 



for all t G [0, r ]. We show now that also the matrix elements U a p(T\ r) and U a ^(T; r) are 
bounded. Indeed, consider 



B(T;r) = U(T;r)U^T;r) . 

We can diagonalise this matrix 

S{T; t)B(T; t)S\T; t) = diag (6 X (T; r), . . . , b M (T; r)) = B diag (T; r) 
S(T;t)S\T;t) = 1. 

The eigenvalues 6 a (T; r) must be positive 

b a {T; r) > , for all r G [0, r ] , (a = 1, . . . , M) . 



Now we consider Tr B(T; r) and find from ( IB. 171) 

9 



<9r 



Tr B(T; r) = Tr {A(T; t)B(T; t) + B(T; t)A\T; t)} 

= Tr {B di , g (T; r) [S(T; r) (A(T; r) + A^T; r)) S+(T; r)] } 

M 

= J2UT;r)A' aa (T;r) , 



a=l 



where 



A' aa (T; r) = [S(T; r) (A(T; r) + A\T; r)) S^T; t 

Suppressing the arguments (T; r) we have 

M 



'cry i 



(9,7=1 
M 



^ (-4/3 7 + A *p) S a pS i 



M 



cry 
A* 



\(9,7=1 / (9,7=1 

A/ 

= e kw+^r 

(9,7=1 

2 



< 4M 2 



2 / ^11 

TO 



(B.25) 

(B.26) 
(B.27) 

(B.28) 



(B.29) 



(B.30) 



(B.31) 



(B.32) 
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Thus we find 



and from (IB~29l) 



\A' aa \<2M-± (B.33) 

T~0 



8 r M 

-Ir5(T;r)<2M^V& tt (T;r) 

dr r ° S (B.34) 

= 2M— Tr 5 (T;r) . 

7~0 



From (1B.34I) and the initial condition (1B.19|) we see that 



Tr Uj{T;t)U\T;t)\ < Mexp(2MC n ) (B.35) 

for all r e [0, To]. This shows that all matrix elements of U(T;t) stay bounded and, in 
particular, have no terms growing exponentially with T. Using (IB. 231) it is easy to show 
that the same holds for all matrix elements of Z7~ 1 (T'; r). 

Consider now an arbitrary vector function r — > £(r) as in (IBID and set 



Cr(r) = ((l-LW))- 1 c) (r) 
= C(r) + (l^Ct) (t) • 



(B.36) 



We have then 



Ct(0) = C(0) , (B.37) 
^Ct(t) = ^ + A(T; r)Cr(r) • (B.38) 



The solution is 



Ct(t) = C(r) + ?7(T; r) f dr' U~ l {T] r')A(T; t')C(t') . (B.39) 

./o 

Since the matrix elements of A(T; r), U(T; r) and f/ _1 (T; r) are all bounded, see (IB.24[) ff. . 
we see that with the norm as in ( IB. 31) we get 

ll&ll < c-IICII • (B.40) 

Here C is a constant. That is, we have from (IB. 361) 

||(1-L (1 ' 1) ) _1 C||<C'||C|| . (B.41) 
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The next term to consider is the following operator occurring in (IB. 141) 

oo 

L (l,2) ^ _ L (2,2))-1 L (2,1) =^L^ (L^Y 

r=0 

OO 

= 5> (r) , 

r=0 

where 

K (r) =L d,2) ( L (2,2)^ L (2,1) _ 



Let £(r) be as in (IB. II) . We get then 

(^ 0) t)jT) = (L^L^()jT) 



N M 



EE'*- 



f3=M+l *e=l 



ri exp 



T T 

i—<p a (n) - i—fein) 

TO T 



x a a p(n) exp [-i7 aa (ri) + vypp{n)] 



n 



x i / dr 2 exp 
'o 



T T 
i—<pp{r 2 ) - i—(Px(T 2 ) 

TO T 



x a /3 ^(r 2 ) exp [-ijpp{r 2 ) + i7**r(T 2 )] G<( r 2) 
(a = l,...,M) . 

Using (14.181) . (IB.lOj) and (IB.llj) we find the following estimates 



N M 



(K^()jT)\<±C 12 C 21 \\C\\AT^ n E/Vexp 

T ° P=M+1 *e=l ^° 



(B.42) 



(B.43) 



(B.44) 



T 



x / dr 2 exp 
'o 



T 
2t~o 



{pp(t 2 ) - pi(r 2 )) 
-i 2 



[r(/3,T 2 )-r(i,r a )] 



(B.45) 



< M(N - M)C 12 C 21 AT^ T 
\\K<®C\\ < M(N - M)C 12 C 21 Ar^|||C|| . 



(B.46) 



In a similar way we get for 1 < a < M 

(K (1 k)jT) = (L^L^L^()jT) 



N N M 



E ED * 



P=M+l 7=M+1 k=\ 



r i exp 



T T 

i— fain) - i— (ppin) 

T~0 T~0 



x a a( g(n) exp [-i7 aa (ri) + i7/9/?fa)] 



x i / dr 2 exp 
Jo 



T T 

i— ¥»/s(Ta) - i— ^ 7 fa) 
r r 



x 0/3 7 fa) (1 - 5/j 7 ) exp [-i7^(r 2 ) + 1777 fa)] 



x i / dr 3 exp 
Jo 



' T T 

i— <P 7 fa) - i— <P*fa) 

To T 



x a 7 ^(r 3 ) exp [-i7 77 fa) + i7^fa)] C*fa) 



AT N M 



(K m ()Jr)\< E E E/ *i«p 



/3=M+1 7=M+1 x=l 



T 
2^ 



fafa) -psfa)) 



n 

x / c?r 2 exp 
'0 

x / dr 3 exp 
Jo 



T 

5— (P/jfa) -P 7 fa)) 

T 

7T- (P 7 fa) "Pi fa)) 
ZTq 



C*12 

r(/?,r 2 )-r(i,r 2 ) 

r oAr min 
r(7,r 3 )-r(i,r 3 ) 
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r Ar r 



a 



21 



Thus we get 



< M(iV - M)C 12 C 21 Ar m J n |(iV - M)C 22 AT^ 



\\KW(\\ < M(N - M)C l2 C 2l AT^(N - M)C 22 AY^ 



It is easy to see that this can be generalised to 

-1 2 



\\K^(\\ < M(N - M)C 12 C 21 AT; 



(iV - M)C 22 Ar, 



-1 1 

lin 



(r = 0,l,2,...) . 

Now we can go back to (IB. 421) . For £(r) as in (IB.lj) we get for large enough T 



|| L (l,2) (1 _ L (2,2) r l L (2,l) C || = \\J2K^( 
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Here we use that for large enough T, ^2^L Q H-K^CII is bounded by a convergent geometric 
series due to (IB. 501) . The analogous argument was already used in (1 A. 71) . 
In an analogous way we can estimate 



(1 _ L (2,2) ) -l L (2,l) C = ^ (L (2,2) ) r L (2,l) C 

where r — > £(r) is as in llB.lj) . We find 

|| (L (2,2) ) r L (2,i) C || <C 2l {TAT mu r^N - M^AT-^Y 

where 



(B.52) 



C 2 \ = 2MC 2 \ max | exp {vy aa {r)) \ . 

M+l<a<N 
re[0,r ] 

With the argument of the geometric series we get, therefore, for large enough T 



(B.53) 



(B.54) 



;i-L( 2 - 2 ))-W)(|<o(i) 



(B.55) 



Consider next the operator L flB.141) . From (1B.51|) and flB.41j) we get, always for large 
enough T, 

iilcii = \\l^ 2 \i - L^y'L^ii - L^rwi 

Iterating this and using again the argument based on the geometric series leads to 



^B.56) 



£i'c|<o(i) 

k=i v 7 



(B.57) 



The next term to consider in ( IB. 13[) is (1 — L^ 2,2 -*) . We analyse this term in a way 
similar to (1 - L^ 1 ))" 1 above, see (jBT5]) - p^5|) . Let 



(B.58) 



\ VcN J 



be a constant vector and set as in 115. lip 
fj a (T; t) = exp 

ae{M+l,...,N} . 



T T 

i—<fa(r) + — Pi(r) + i7 aQ (r) 

To 2Tq 



(B.59) 



We define a matrix function V(T; r) by 



TV 



V a p(T; t) VcP = ((1 - L^)-^) q (T; r) 

I3=M+1 



This matrix satisfies 



V(T;r) = [y aP (T;r)) , 
a,/3 e {M + 1, . . . , N} . 

V{T; 0) = 1 N - M , 
^-V(T;t)=A(T;t)V(T;t) , 



where we get with ( 14. 3ft - (14.61) 

A(T-t)= [A a(3 (T-T 
A al 3(T;r) = 

a,@e {M + 1,...,N} 
In analogy to (1B.25j) ff. we consider 



i-E(a,r) + ^r(l,r) 

To 2T 



5 a /3 + iOa/3( 



B(T;r) = V(T;r)V^T;r) 



and get 



|- Tr £(T; r) = Tr {b(T; r) [l(T; r) + i+(T; r) 



< 



r AT , 2c 

'- 1J - mm ~ 

To T . 



Tr5(T; 



where we used ( 15. 2 p and set 



r 



2 TG[0,T ] 



TV 



1/2 



max |a Q/ 3(r) - a* 0a (r) 



From (IB.66P we get easily 



T 



TiB(T;T)<(N-M)e 2c ex V ^--AT min Tj . 
\Vap(T- r)| < v/iV^M e c exp (-I-AT^t^ 
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for all a, P E {M + 1, . . . , N}. With the Cauchy-Schwarz inequality we get from (1B.68j) also 

-i 1/2 

E i^-)i 2 )( E ^ 



N 



E 1^(^)1 < 

a,/3=M+l 



va,/3=M+l 



1/2 



= (N-M) p(T;r) 
< (JV - Mf/Vexp ^-I-AI^r) . 



From (IB. 701) we find 

m-L^)- l v)Jr)\ 



N 



E V aP {T)T)r} c t 
T 



< C 22 exp ( -— Ar min r , N , /; N . 



a E {M + 1, . . . , N} 



where 



C 22 = {N- Mf 2 e c . 
With the result (IB. 701) we can estimate the term 

£(1,2) (1 _ L (2,2))~l£(0) = L (l,2)^(0) 



in (1B7L3]) . where £ (0) is given in Q and £ (0) in f loTTTD . We get from flB~6l) 



(LMVt&>) a (r) = i rfr'exp 



o 



T 



T 



i—ip a (j') - — pi(r') - i7 aQ (r') 



to 



2r 



A' 



x E MWioef, 

/3, 7 =Af+l 



a G {1,...,M} ; 
|(L(^V-^) a (r)| <^||^|| jT^ E 

where we set 



A? 



/3, 7 =M+1 



C' 12 = max r |a a/3 (r) exp (-i7 Qa (r)) | . 

l<a<M 
M+l<p<N 
re[0,T ] 



(B.70) 



(B.71) 



(B.72) 



(B.73) 



(B.74) 



(B.75) 



(B.76) 
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With (lB~70l we obtain finally 



|| L (i 1 2) (1 _ jL (2,2) ) -i|(0)|| = || L d,2)^(0)|| < ±AT^ m C[ 2 (N - M) 3 / 2 e ( 



(B.77) 



Now we have collected all tools needed to discuss the solution for x(T] t) given in (I5.17P 
and (IB. 13|) for large T. We have from (1B.13|) 



X = (1 - [x (0) + L<M (1 - L (2 ' 2) ) f (0) 

oo 

+ (i - z^ 1 ' 1 )) -1 lk [x i0) + l m (i - L^y 1 i 

k=l 

From (IB. 77ft we have 

l| L (i, 2) (i - lM)- 1 = o \\t®\\, 

With (JR57D and (lB~iT|) we find 



(o) 



fc=l 



(i-i (1 -T 1 E^ (1 ' 2) (i-^ a,2, )" 1 e (0) || = o(^) ne (0) n . 
fc=i ^ ' 

Inserting all this in (1B.78|) gives 

With the definition of the matrix U(T; r), ( IB. 161) . we get explicitly 
X (T; r) = U(T; r) X (0) + O (^) ||x (0) II + O (± 



(B.78) 



(B.79) 



(B.80) 
(B.81) 



(B.82) 



(B.83) 



[1] M. DeKieviet, D. Dubbers, C. Schmidt, D. Scholz, and U. Spinola, Phys. Rev. Lett. 75, 1919 
(1995). 

[2] Y. B. Zel'dovich, Sov. Phys. JETP 9, 682 (1959). 

[3] M. A. Bouchiat and C. Bouchiat, J. Phys. (Paris) 35, 899 (1974). 

[4] M. Bouchiat and C. Bouchiat, Rep. Prog. Phys. 60, 1351 (1997). 

[5] J. Guena, M. Lintz, and M. Bouchiat, Mod. Phys. Lett. A 20, 375 (2005). 



30 



[6] I. B. Khriplovich, Parity Nonconservation in Atomic Phenomena (Gordon and Breach, 

Philadelphia etc., 1991). 
[7] S. C. Bennett and C. E. Wieman, Phys. Rev. Lett. 82, 2484 (1999). 

[8] P. A. Vetter, D. M. Meekhof, P. K. Majumder, S. K. Lamoreaux, and E. N. Fortson, Phys. 

Rev. Lett. 74, 2658 (1995). 
[9] M. V. Berry, Proc. R. Soc. Lond. A392, 45 (1984). 
[10] A. Shapere and F. Wilczek, eds., Geometric Phases in Physics, vol. 5 of Advanced Series in 

Mathematical Physics (World Scientific, Singapore, 1989). 
[11] T. Bergmann, T. Gasenzer, and O. Nachtmann, Metastable states, the adiabatic theorem and 

parity violating geometric phases II, HD-THEP-07-09, physics/0703276 (2007). 
[12] G. W. Botz, D. Bruss, and O. Nachtmann, Ann. Phys. (NY) 240, 107 (1995), hep-ph/9406222. 
[13] A. Messiah, Quantum Mechanics, vol. 2 (North Holland, Amsterdam, 1970), 6th ed. 



